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Abstract 

During the spin-down phase of the life of a higher-dimensional black hole, the emission 
_ of particles on the brane exhibits a strong angular variation with respect to the rotation 

axis of the black hole. It has been suggested that this angular variation is the observable 
that could disentangle the dependence of the radiation spectra on the number of extra 
dimensions and angular momentum of the black hole. Working in the low-energy regime, 
we have employed analytical formulae for the greybody factors, angular eigenvalues and 
eigenfunctions of fermions and gauge bosons, and studied the characteristics of the cor- 
responding angular profiles of emission spectra in terms of only a few dominant partial 
modes. We have confirmed that, in the low-energy channel, the emitted gauge bosons 
become aligned to the rotation axis of the produced black hole while fermions form an 
angle with the rotation axis whose exact value depends on the angular-momentum of the 
black hole. In the case of scalar fields, we demonstrated the existence of a "spherically- 
symmetric zone" that is followed by the concentration of the emission on the equatorial 
plane, again in total agreement with the exact numerical results. 



1 Introduction 



Under the assumption that a low-energy scale for gravity exists in the context of a higher- 
dimensional fundamental theory [Tj, the possibility of observing in the near future quantum- 
gravity effects has excited a lot of interest among high-energy physicists, both theorists and 
experimentalists. The main reason for that is the fact that, if M* - the fundamental gravity scale 

- is as low as a few TeV, then these effects could be observed during trans-Planckian particle 
collisions at current, ground-based accelerators [2]. One such strong- gravity effect could be the 
creation of higher- dimensional miniature black holes during the collision of ordinary Standard- 
Model particles localised on our brane - a (3+l)-dimensional hypersurface embedded in the 
(4 + n)-dimensional spacetime, the bulk. Due to their small size, these black holes will have a 
high temperature and will evaporate very quickly via Hawking radiation [3], i.e. the emission 
of ordinary particles with a thermal spectrum [H El |6]. 

The emission of Hawking radiation is anticipated to take place during the two intermediate 
phases in the life of the black hole, the spin-down and the Schwarzschild phase. It is expected to 
be the main observable signal not only of the creation of these miniature black holes but of the 
existence of the extra spacelike dimensions themselves in the absence of which the creation of 
the former would not be possible. As a result, the study of the emission of Hawking radiation 
by higher-dimensional black holes has been intense during the last ten years. In the early 
days, the Schwarzschild phase - the spherically-symmetric phase in the life of the black hole 
arising presumably after the shedding of its angular momentum - was considered to be the 
longest and thus the most important. It was also the one with the simplest metric tensor 
describing the spacetime around it, and therefore the first one to be exhaustively studied both 
analytically [TJ [8] and numerically [9] [10]. The results derived showed a strong dependence 
of the emission rates of all types of Standard Model particles on the brane on the number of 
spacelike dimensions existing transversely to the brane 0. 

One was thus led to hope that by detecting the emission of Hawking radiation could not 
only shed light on aspects arising from the interplay between classical gravity and quantum 
physics but also give a quantitative answer to a century-old fundamental question, that of the 
dimensionality of spacetime. Nevertheless, the Schwarzschild phase is preceded by the axially- 
symmetric spin-down phase. The gravitational background around a simply-rotating black hole 

- one of the very few cases where the equations of motion of the propagating particles can be 
decoupled and solved - depends also on the angular-momentum parameter a of the black hole. 
According to the results existing in the literature [IlEHEIllIIllIHlIISI^ 

this dependence is carried over in the form of the radiation emission spectra and is, in fact, 
found to be similar to the effect that the number of additional spacelike dimensions n has on 
them. To complicate things more, simulations of black hole events [27J [28] have revealed that 
the spin-down phase is not a short-lived one, as previously thought, and that the rotation of 
the black hole remains significant for most of its lifetime. 

The fact that the dependence of the radiation spectra on the number of extra dimensions n 
for all types of particles is entangled with the dependence on the angular-momentum parameter 
a means that measuring both of these parameters is extremely difficult. The only way out was 

Variants of the spherically-symmetric Schwarzschild phase, where a cosmological constant llh or the higher- 
curvature Gauss-Bonnet term |12j were introduced, were also studied with the spectrum exhibiting a dependence 
also on parameters related to these terms. In addition, the Schwarzschild phase of quantum-corrected black 
holes has been studied in [T5] . 
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to employ another observable that would strongly depend on only one of these two parameters 
while being insensitive to the other. Upon determination of that particular parameter, the 
second could then be determined from the radiation spectra. One characteristic feature of the 
emission spectra coming from the spin-down phase is the non-isotropic emission, in contrast to 
the one coming from the Schwarzschild phase where the emitted particles are evenly distributed 
over a 4tc solid angle. It has therefore been suggested [291 [30] that this non-isotropy can serve 
as the additional observable necessary to disentangle the n and a-dependence of the spectra. 
Indeed, it was demonstrated [30] that the angular profile of the emitted radiation depends 
extremely weakly on the number of additional dimensions n while it may provide valuable 
information on the angular momentum of the black hole (see, for example, [31]). 

More specifically, under the combined effect of the centrifugal force exerted on the emitted 
particles and the spin-rotation coupling for particles with non-zero spin (an analytical explana- 
tion of the latter is given in [32]), the orientation of the emitted radiation depends strongly on 
the energy channel in which the particles are emitted and on how fast the black hole rotates. If 
we look specifically at the low-energy channel, then we observe that gauge bosons and fermions 
have a distinctly different behaviour: the emitted gauge bosons remain aligned to the rotation 
axis of the black hole independently of the angular-momentum parameter; fermions, on the 
other hand, form an angle with the rotation axis whose value strongly depends on the value of 
a. As a result, the orientation of gauge bosons can serve as a good indicator of the rotation 
axis of the black hole [30] and the orientation of fermions can then provide a measurement of 
the value of the angular momentum of the black hole [29], [30] . 

The aforementioned results presented in [291 130] were derived by means of a very compli- 
cated and time-consuming process that involved the numerical integration of both the radial 
and angular part of the equation of motion of each emitted particle as well as additional chal- 
lenges such as the numerical calculation of the angular eigenvalue itself, which does not exist 
in closed form for a rotating background, and the summation of a very large number of partial 
modes. The purpose of this work is to provide an alternative way of deriving the angular profile 
of the emitted radiation without resorting to complicated numerical calculations. This is facil- 
itated by the fact that all valuable information that may be derived from the angular spectra 
is restricted in the low-energy regime where the radial equations for all types of particles have 
been analytically solved [T8| IT9] . In addition, analytical formulae, in the form of power series, 
for the angular eigenfunction and eigenvalue exist in the literature. By combining all the above 
in a constructive way, we investigate which contributions are the dominant ones, that predomi- 
nantly determine the angular profile of the emitted radiation. In this way, we formulate simple 
constraints involving a finite number of terms and partial modes that successfully reproduce 
all the features of the anisotropic emission, namely the value of the angle where the emission 
becomes maximum and the corresponding value of the energy emission rate. 

The structure of this paper is as follows. In section [2j we present the theoretical framework 
with the field equations that need to be solved and the corresponding energy emission rates for 
a general spin-s field. In section [21 we present the analytical formulae for the greybody factors, 
angular eigenfunctions and eigenvalues that will be our tools for the analytical investigation of 
the angular profile of the emitted radiation. In section HI we consider separately the cases of 
fermions, gauge bosons and, for completeness, scalar fields too, emitted by a higher-dimensional 
simply-rotating black hole on the brane: in each case, we determine the dominant modes, 
formulate simple extremization constraints with respect to the angle of emission 9, and derive 
their angular distribution on the brane. Finally, in section El we summarise our results and 
present our conclusions. 



2 Theoretical framework 



The most generic type of a black hole in a higher-dimensional spacetime is the one that rotates 
around one or more axes. The gravitational field around such a black hole is described by 
the Myers- Perry solution [33] . However, it is only for particular configurations of the angular- 
momentum components that the equation of motion of a particle propagating in the higher- 
dimensional spacetime can be decoupled into an angular and a radial part. The case of a simply- 
rotating black hole, where the black hole possesses only one angular-momentum component that 
lies on a plane parallel to our brane, corresponds to one of these configurations and the one that 
has been mostly considered in the literature. This choice is also justified by the assumption 
that the black hole, if created by the collision of two brane-localised particles, will acquire an 
angular momentum component along the (3+l)-dimensional part of the full manifold. 

In this work, we will also focus on the case of a simply-rotating black hole. In addition, 
we will study effects that take place strictly on our brane, namely the emission of Hawking 
radiation by the higher- dimensional, rotating black hole in the form of non-zero-spin Standard- 
Model fields. The line-element of the brane background in which these particles propagate is 
given by the expression pE] 

ds 2 = - ( 1 - -JL-) - 2a " ^ 9 dt dip + (r 2 + a 2 + aV Sin ' 6 ) sin 2 9 V 

+ ^dr 2 + Ed# 2 , (1) 

where 

A = r 2 + a 2 --^, £ = r 2 + a 2 cos 2 #. (2) 

The mass Mbh of the black hole and its angular momentum J are related to the \i and a 
parameters, respectively, through the relations 

(n + 2)A n+2 2 

M BH = t—^ H, J = — —;M BH a , (3) 

IbnG n + 2 

where A n+2 = 2ii ( - n+3 ^ 2 /T[(n + 3)/2] is the area of an (n + 2)-dimensional unit sphere, and G 
is the (4 + n)-dimensional Newton's constant. The horizon radius follows from the equation 
A(r) = 0: for n > 1, it may be shown that there is only one real, positive root, which may be 
implicitly written as r^ +1 — + a 2 ), where a* is defined as a* = a/r^. 

The derivation of the field equations that the brane-localized Standard- Model fields satisfy 
in the above background follows the analysis performed originally by Teukolsky in 4 dimensions 
[34J . The method demands the use of the Newman- Penrose formalism and results in a 'master' 
partial differential equation that scalars, fermions and gauge bosons obey on the brane. If we 
use a factorized ansatz for the field perturbation ^ of the form 

* h (t, r, 6,p) = J2 ha A h R A (r) h S A (6) e** , (4) 

A 

the aforementioned 'master' equation separates, in the background of Eq. (CQ), into two decou- 
pled ordinary differential equations, a radial 
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l S A (x) = 0. (6) 



In the above, h is the spin-weight, h = (—\s\,+\s\), of the given field that distinguishes its 
radiative components, and A = {Imu} denotes the set of 'quantum numbers' of each mode. 



We have also defined the quantities K 



+ a )tu — am and x = cos#. Finally, hA\ is 



the eigenvalue of the spin- weighted spheroidal harmonics hSh(x) - as we will shortly comment, 
the value of this constant does not exist in closed form. This quantity also determines the 
separation constant between the radial and angular equations with ^Aa = hA\ — 2mau + a 2 u 2 . 

The above set of equations has been used in the literature in order to study the emission of 
Hawking radiation, in the form of an arbitrary spin-s field, from a higher-dimensional, simply- 
rotating black hole on the brane [HI [151 E] ■ The resulting differential energy emission rate per 
unit time, energy and angle of emission is given by the expression [TU [T5] 
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The radiation spectrum of the black hole resembles those of a black body with a temperature 

(n + l) + (n-l)al 
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At the same time, however, the spectrum is significantly modified compared to the black-body 
one: in the exponent, the combination £j = u — am/ (a 2 + r 2 ), includes the effect of the rotation 
of the black hole; also, the quantity Ta, the transmission probability (or, greybody factor), 
determines the number of particles that eventually overcome the gravitational barrier of the 
black hole and reach asymptotic infinity. If Eq. ([7]) is integrated over all angles of emission 6, 
we obtain the power rate in terms of unit time and energy 
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(9) 



The derivation of the integrated-over-all-angles power spectra, for all species of brane- 
localised fields - scalars, fermions and gauge bosons, was performed both analytically [HI [19] 
and numerically [TU [TSJ [T7] . According to these results, the energy emission rate - as well as the 
particle and angular-momentum emission rates - are significantly enhanced as both the number 
of additional, spacelike dimensions and the angular-momentum of the black hole increase. The 
enhancement factor was of order (9(100) when n varied between 1 and 6, and of order O(10) 
as a* increased from zero towards its maximum value a™ ax — (n + 2)/2. 

In contrast to the case of the spherically-symmetric Schwarzschild phase, the emission 
of particles during the rotating phase of the life of the black hole is not isotropic. The axis 
of rotation introduces a preferred direction in space and the emitted radiation exhibits an 
angular variation as 9 ranges from to n. It was found [TU [TSJ [T7] that a centrifugal force is 
exerted on all species of particles, that becomes stronger as either u or a increases and forces 
the particles to be emitted along the equatorial plane (9 = n/2). In addition, for particles 
with non-vanishing spin, an additional force, sourced by the spin-rotation coupling, aligns the 
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emitted particles parallel or antiparallel to the rotation axis of the black hole - this effect is 
more dominant the smaller the energy and larger the spin of the particle is. If the form (J7J) of 
the power spectrum is used where both helicities appear, the spectrum is symmetric over the 
two hemispheres, 9 G [0, 7r/2] and [tt/2, tt]. If a modified form, in which only one of the helicities 
appear each time, is used instead, then the angular profile is asymmetric with particles with 
positive helicity (corresponding to -\ a \S%) being emitted in the upper hemisphere and particles 
of negative helicity (corresponding to + | s |Sa) being emitted in the lower one. This angular 
variation in the Hawking radiation spectra is considered to be one of the main observable 
effects on the brane of a higher-dimensional, rotating, decaying black hole. 

One would ideally like to deduce the values of both spacetime parameters, n and a, from 
the predicted forms of the Hawking radiation spectra. However, the fact that both parameters 
affect the integrated-over-all-angles spectra in a similar way impose a great obstacle. The 
resolution of this problem would demand the existence of an observable that depends strongly 
on only one of the two parameters while being (almost) insensitive to the other. That observable 
was shown [29, 30j to be the angular variation of the spectra discussed above. Particularly, in 
the low-energy channel, the alignment of the gauge bosons along the rotation axis can reveal 
the orientation of the angular-momentum of the black hole. Then, it was demonstrated that 
the angle of emission of fermions, in the same energy channel, is very sensitive to the value of 
the angular-momentum of the black hole: the larger the a parameter is, the larger the value of 
9, around which the emission is peaked, becomes. Remarkably, the behaviour of gauge bosons 
and fermions alike remains unaltered as the dimensionality of spacetime changes. 

3 Analytical forms of the radial and angular functions 

The results on the angular profile of the emitted fields with non-zero spin on the brane, discussed 
above, were derived by numerically integrating both the radial (JSJ) and the angular (jSj) equation: 
the latter in order to find the angular eigenvalue hAa and eigenfunction ^Sa, and the former 
in order to determine the greybody factor Ta through the radial function hR\. The numerical 
manipulation of the radial and angular differential equations is necessary for the derivation of 
the exact solutions for hR^ and hS\, respectively, and subsequently of the complete Hawking 
radiation spectra. However, when it comes to the spectra of gauge bosons and fermions revealing 
information about the orientation of axis and value of the angular momentum of the black hole, 
the range of interest is the low-energy one. Thus, in what follows we will focus on the low-energy 
channel, and attempt to derive analytically information about the angular profile of non-zero- 
spin fields emitted on the brane. To this end, we will henceforth ignore the single-component 
scalar fields and concentrate our study on brane- localised fields with spin 1/2 and 1. 

Under the assumption of low-energy of the emitted field and low-angular-momentum of the 
black hole, the radial equation (JSJ) was analytically solved in [HI [19] for all species of particles. 
A well-known approximation method was used in which the radial equation was solved first near 
the horizon, then at asymptotic infinity, and the two were finally matched at an intermediate 
regime to construct the complete solution for hR\- The transmission probability Ta for fermions 
was defined as the ratio of the flux of particles at the black-hole horizon over the one at infinity, 
with the flux being determined through the conserved particle current. For gauge bosons, 
where no conserved particle current exists, a radial function redefinition and a simultaneous 
change of the radial coordinate conveniently change the corresponding gravitational potential 
to a short-range one - then, the amplitudes of the outgoing and incoming plane waves at infinity 

K 



can easily determine the transmission probability. For fermions and gauge bosons, Ta comes 
out to have the form 
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In the above, the quantity Z, defined by 

Z= y/(2\s\ -l) 2 + 4( h A A + 2\s\ + a 2 co< 



(13) 
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appears in the solution of the radial equation in the asymptotic infinity that is expressed in 
terms of the Kummer functions M and U. Similarly, the coefficients (a, b, c), given by 



a = a + (3 + — 1 



b = a + (3 . 



2a 



(15) 



are the coefficients of the hypergeometric function F in terms of which the solution of the radial 
equation is written near the black-hole horizon. Finally, the following definitions hold 
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supplemented by the following ones: A* = n + 1 + (n — l)a 2 , = K/r h and h^A = h^A + 2|s|. 

For scalar fields, the transmission probability is again defined from the amplitudes of the 
outgoing and ingoing spherical waves at infinity [18] 

B — i 
B + i 

where B now is given by the expression 
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(19) 



We note that the angular eigenvalue hA\ makes its appearance in the above analytic 
results both in Eq. (THE)) and Eq. (II 7p . As already mentioned in the previous section, in the 
case of a rotating black hole, this quantity does not exist in closed form. For arbitrary large 
values of the energy of the emitted particle and angular momentum of the black hole, its value 
can be determined only via numerical means - that was the method applied in [T4"l [15] where 
the complete spectra for scalars, fermions and gauge bosons were derived. However, for low u 
and low a, the angular eigenvalue of the spin-weighted spheroidal harmonics can be expressed 
as a power series with respect to au [35j EHl E7J ES [39] 
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By using the above power-series form for the angular eigenvalue and keeping terms up to 
fourth order, the analytically derived formulae for the transmission probabilities (FIU]) and (FiT|) 
for fermions and gauge bosons - as well as the one for scalar fields - were shown in [181 EE] to 
be in excellent agreement with the exact numerical ones derived in [HI H3) . The power-series 
expansion of the angular eigenvalue is quite cumbersome and, up to the sixth order, can be 
found in [23 ESI E3 ESI EH]- It is worth giving here, some particularly simple formulae we have 
derived, for the needs of our analysis, for the eigenvalues of fermions and gauge bosons up to 
second order, namely 
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max (|m|, \s\) 

In the above, we have given the values of the angular eigenvalues for the positive helicities 
h = 1/2 and h = 1, respectively. The angular eigenvalues exhibit a well-known symmetry 
[40l [39] according to which, if \ s \Aa is the eigenvalue for the positive-helicity component of 
a given field, then the one for the negative helicity _| s |^4a readily follows from the relation 



\s\A\ = \ s \ 



A\ + 2\s\. For completeness, we add here a similar formula for the angular eigenvalue 



of scalar fields that first appeared in [6]: 
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(24) 



Let us now turn to the angular equation ()6|). Leaver [40J found an analytic solution for 
the angular eigenfunction hS\(x) that may be expressed as a series of the following form 



^(x) = e^ x (1 + x) k - (1 - x) fc + a p (1 + xy 



(25) 



p=0 
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where x = cos9 and k± = \m ± /i|/2. The expansion coefficients a p can be found through a 
three-term recursion relation 

a a x + /3 a = , (26) 
a p a p+1 + (3 P a p + 7 P a p _i = , 0=1,2,...) (27) 
In the above, the coefficients (a p , (3 P , j p ) are in turn determined by the relations 

a p = -2(p+l)(p + 2A;_ + l), 

A> = p(p - 1) + 2p(fc_ + k + + 1 - 2aw) - [aV + + 1) + h A A ] 

- [2aco(2k_ + h + l) ~ (k_ + k + )(k_ + k + + l)], (28) 

7 p = 2au(p + k_ + k + + h) . 

The above analytic form determines the angular eigenfunction up to a constant that can be fixed 
by imposing the normalization condition J \hS\(x)\ 2 dx = 1. According to [39], an excellent 
approximation to the exact solution is obtained by keeping ~ 10 terms in the expansion of (T25]) . 



4 Analytical description of the angular profile 

In this section, we will attempt to study the angular profile of the emitted Hawking radiation 
on the brane by employing semi-analytic techniques. Our starting point will be Eq. ([7]) that 
determines the angular profile of the emitted radiation as a function of x = cos 9. By using the 
analytical formulae presented in the previous section, we will compute the value of the angle 9 
where the emission of particles becomes maximum. Since the emission of positive and negative 
helicity components is symmetric under the change 9 — > n — 9 [30], in what follows we consider 
only the emission of positive helicity components, h > 0. 

In Eq. ([7]), the dependence on the angle 9 is restricted in the angular eigenfunction h SA(x). 
One may then naively try to extremize this equation to find a constraint that will determine 
the desired value of 9 max , defined as the value of the angle where the differential rate of emission 
takes its maximum value. We then obtain 

dx [dxdtdoj} 47r exp(a}/T^) ± 1 A ( h A dx ) ^ ^ 

By employing the analytical expression (|25|) for the angular eigenfunction hS\{x) and evaluating 
the derivative, we obtain the following constraint 

J2 h W A (l + x) 2k -(l-x) 2k + ^a p (l + x)^ ffl! (l + xW flW + T -^- T -^ =0. (30) 

l,m P=0 q=0 V 7 

In the above, we have defined the "weight factor" hW\ as 

* W *= f? OUT U1 Ta " (31) 

exp[(w — mil) / 1 h) ± 1 

The analytical evaluation of the constraint ( 130]) in full is not possible. As mentioned 
above, the sum over p (and q), originating from the analytic form of the angular eigenfunction, 



may be truncated at a finite value, but care must be taken so that the truncated series remains 
close to the exact solution and the value of # max is not affected. The constraint contains two 
additional sums: one with respect to /, the total angular- momentum number ranging from \s\ 
to oo, and one over m, the azimuthal angular-momentum number that takes values in the range 
[—/,+/]. None of these sums can be discarded: all of the quantities involved, the coefficients 
k±, a p (and a q ), as well as the weight factor hWa, depend on both angular-momentum numbers 
in a non-trivial manner. It is, therefore, the combined contributions of all, in principle, partial 
modes that determines the angular profile of the emitted radiation. Finally, these contributions 
do not enter on an equal footing: each mode carries a weight factor hWh ~ defined in Eq. (I3TI) 
in terms of the 'thermal/statistics' function and the greybody factor Ta - that determines the 
magnitude of its contribution to the angular profile. 

In what follows, we will attempt to shed light to the important contributions to Eq. (!30l) 
that determine the value and location, in terms of the angle 9, of the maximum emission rate 
for fermions and gauge bosons. As the interesting phenomena take place in the low-energy 
regime, we will use purely analytic expressions for all quantities involved, namely the angular 
eigenvalue, the angular eigenfunctions and the greybody factor. Having been established in the 
literature [30] that the orientation of the emission of fermions and gauge bosons is not affected 
by the value of the number of extra dimensions introduced in the model, we will keep fixed the 
value of n and, henceforth, set n = 2. 

4.1 Emission of Fermions 

We will start with the most phenomenologically interesting case, the emission of fermions. Our 
strategy will be the following: by using the most complete analytic forms, we will investigate 
when a particular contribution to the angular profile becomes so small that is irrelevant and 
can thus be ignored. We will therefore use the power series expansion (1201 for the angular 
eigenvalue up to fourth order in (acu), the analytic form of the angular eigenfunction given in 
f|25l) by keeping terms up to p = 10, and, at the beginning, allow the angular-momentum 
numbers (7, m) to vary over their full range. 

In Figs. H](ab), we depict the differential emission rate ([7]) per unit time, unit frequency 
and angle of emission in terms of cos 8, for the case = 0.5 and a* = 0.5 (left plot) and 
a* = 1.5 (right plot). The different curves correspond to the derived spectrum where modes 
up to a certain value of I (and all values of m in the range [—1, +1] ) have been summed up: 
the lower (blue) curve includes only the 1 = 1/2 modes, the next (green) one modes up to 
/ = 3/2, the subsequent (red) one modes up to I = 5/2 and the last (orange) one modes up to 
/ = 7/2. We observe that the 1 = 7/2 curve is not even visible as it is completely covered by 
the 1 = 5/2 one - the same happens for all higher modes. As a matter of fact, the difference 
between the 1 = 5/2 and 1 = 3/2 curves is also quite small: for the maximum value of the 
angular momentum considered, a* = 1.5, the difference in the value of the emission rate at its 
maximum and of # max is of the order of only 1%; for smaller values of a*, the errors reduce even 
more: for a* = 0.5, the difference in the value of the emission rate at its maximum drops at 
the level of 0.08% while 6* max is not affected at all. We may thus conclude, that the sum over 
I in (1301) can be safely truncated at I = 3/2. The reason for this significant truncation is the 

2 We have confirmed that, by keeping terms up to p = 10 in this expansion, the derived values of the angular 
eigenfunction agree extremely well with the exact numerical ones - as a consistency check, we have successfully 
reproduced the plots of the angular eigenfunctions appearing in [351 US] . 
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Figure 1: The differential energy emission rate (JTj) in terms of cos#, for n = 2, cu* = 0.5, and 
(a) a* = 0.5 (left plot) and (b) a* = 1.5 (right plot). The different curves correspond to the 
emission rate when partial modes up to Z = 1/2, Z = 3/2, Z = 5/2 and Z = 7/2 (from bottom to 
top) have been summed up. 

weight factor hWx- although the thermal/statistics factor gives a boost to modes with large 
and positive m, the significant suppression of the greybody factor Ta in the low-energy regime 
as Z increases ensures that higher modes can be safely ignored. 

As a next step in our study, we investigate whether the sum in the series expansion 
of the eigenfunction can also be truncated. To this end, we have computed the differential 
energy emission rate (j7]), for = 0.5 and a* = 1, by keeping modes up to Z = 5/2 for extra 
safety, and gradually increasing the maximum value of the sum index p. The behaviour of the 
corresponding results for the emission rate as a function again of cos# is plotted in Fig. |2} 
where the different curves correspond to the maximum value of p kept in the sum, p = 1, 2, 3, 4 
and 10. We observe that the correct value of the emission rate at its maximum is obtained 
fairly soon, when terms only up to p = 2 are included in the sum; the value of 9 max , on the 
other hand, needs one more term in the expansion (p = 3) to acquire its actual value. Our 
results are not in contradiction with [39] where the value of p = 10 was defined as the one that 
accurately reproduces the exact form of the eigenfunction. Indeed, higher terms included in 
the sum up to p = 10 do change the behaviour of the eigenfunction, however, these changes 
are restricted in the area away from the angle of maximum emission, as Fig. [2] clearly shows. 
The value of the angular momentum of the black hole strongly affects the value of p maiX : for 
a* = 1.5, the correct value of 9 is obtained when terms up to p = 4 are included; in contrast, 
for a* = 0.5, no terms higher than p = 1 are needed in the sum. 

Let us comment at this point on the expression of the angular eigenvalue that was used in 
our calculations. As noted above, we initially employed the power series form of Eq. ( T2"Uj) with 
terms up to the fourth order in (au). However, we have found that the expression f T2T|) . with 
terms up to second order only, is more than adequate to lead to accurate results. Although 
including higher-order terms cause, at times, a significant change in the value of the angular 
eigenvalue itself, that change hardly affects any aspects of the angular profile of the emitted 
radiation. For example, for the mode Z = 1/2 and m = —1/2, the difference in the value of 
the eigenvalue, when terms up to second and third order, respectively, are kept, is of the order 
of 10%, the effect in the value of the coefficient (3 P appearing in Eq. f l28|) is only 0.2% which 
leaves the angular profile virtually unchanged. 
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Figure 2: Energy emission rate per unit time, unit frequency and angle of emission in terms 
of cos (9, for n = 2, a>* = 0.5 and a* = 1, and terms in the series expansion of the angular 
eigenfunction summed up to p — 1, p — 2, p — 3, p — 4 and p = 10. 



One may simplify further the analysis by considering more carefully the partial modes 
that dominate the energy emission spectrum. According to the results above, the sum over I can 
be safely truncated at the value I = 3/2, and thus we need to sum over the following six modes: 
(l,m) = [(|, §), (§, — §), (§, f), (§, §), (§, — |), (§, — §)]. However, not all of the above modes 
have the same contribution to the angular variation of the energy spectrum. In Fig. E^a), we 
display the angular eigenfunctions of the four most dominant modes out of the aforementioned 
six, for n = 2, = 0.5 and angular momentum a* = 0.5 (left plot). It is clear that, for 
small values of a*, the two 1 = 1/2 modes dominate over the 1 = 3/2 ones. This dominance 
is further enhanced when the corresponding weight factors are taken into account, with the 
ones for the 1 = 3/2 modes being at least one order of magnitude smaller than the ones for 
the I = 1/2 modes. But even the contribution of the two dominant modes, (~,±~), is not of 
the same magnitude: when the weight factors and the difference in magnitude of the angular 
eigenfunctions are taken into account, the |) mode is found to have at least five times bigger 
contribution than the (|, — |) one. As a result, the angular pattern of the emitted radiation at 
the low-energy channel, for small values of the angular momentum parameter, is predominantly 
defined by the (|, |) mode. Then, the constraint ( I30p takes the simplified form| 
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In the above, we have used that fc_ 



and k + = 1/2 for the mode (|, |), and the superscript 



{1/2} denotes that the set of coefficients a^ 2 ^ for this particular mode should be used here. In 



3 In what follows, we will adopt the value p = 3 as the maximum value of the sum index needed to accurately 
reproduce the behavior of the fermionic eigenfunction around the angle of maximum emission. 
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Figure 3: The fermionic angular eigenfunction i/2<S/ m as a function of cos 9, for n = 2, ou* = 0.5 
and: a) a* = 0.5 (left plot) and (I, m) = [(|, ~), (~, — ~), (§, §), (§, 5)] (from top to bottom), b) 
a* = 1.5 (right plot) and (l,m) = [(|, |), (|, |), (|, ~), (~, — |)] (from top to bottom). 



Appendix lA.lt we list the results for the angular eigenvalue, as this follows from Eq. ( 121 p . the 
values of the (a p ,{3p,j p ) coefficients, according to the definitions (I2"g]) . and finally the relations 
between the first four sum coefficients a p , given by the three-term recursion relations fl26]) -( l2~T]) . 
A simple numerical analysis, then, shows that Eq. ( 133]) does not have any roots in the range 
x G (—1,-1-1) for auj < 0.52, with the global maximum located at x — — 1 and the global 
minimum at x = 1. Therefore, if we fix the energy channel at e.g. = 0.5, the angular 
eigenfunction of the (|, |)-mode does not show any extrema up to a* ~ 1; as a result, the 
energy emission rate takes its maximum value at 9 = tt in accordance with the exact numerical 
results derived in J29J [30] . 

Nevertheless, as a* increases, the (|, |)-mode becomes important - this may be clearly 
seen in Fig. [3]^b). Let us examine the behaviour of this mode on its own. Its extremization 
constraint is given now by 

(3/2) + a (3/2) (1 + x)+ a (3/2) (1 + x)2 + a (3/2) (1 + ^3 

+aj 3/2) + 2af 2) (1 + x) + 3c4 3/2) (1 + x) 2 = , (34) 

where we have used that, for this mode, A;_ = 1/2 and k+ = 1. By making use of the 
relations between the first four a^ 2 ^ coefficients, as these are found again in Appendix lA.l] 
and performing a simple numerical analysis, we arrive at the following results: for au = 0, all 
<2j with i > 1 vanish, and the constraint (134"]) reveals the existence of a sole extremal point at 
x = —1/3; this extremum is a local maximum - as a* increases, the local maximum becomes 
gradually more important and slowly moves to the left, thus competing with the maximum of 
the (|,|)-mode at x — —1 to create a global maximum for the energy emission rate in the 
range (—1, —1/3) with the exact location depending on the value of a*. 

Thus, summarizing the above results, for an arbitrary value of a*, the angular variation of 
the emitted fermions is mainly determined by the contribution of the (~, ~) and (|, |) modes, 
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Figure 4: The relative weight factor W re i = W1/2/W3/2 in terms of the angular-momentum 
parameter a*, for the particular case of = 0.6. 

and thus the constraint fl30|) may take the final form 

w„, £ < m (1 + *? £ 4 1/2) + (<- + - 57^-1, 

p=0 g=0 ^ I A 

3 3 

+(1 - x 2 ) E 4 3/2) (i + *) p E 4 3/2) (i + (<*> + tt^ - t^- ) = • ( 35 ) 

We have also defined the relative "weight factor" W re i = W1/2/W3/2 whose value depends 
strongly on the angular parameter a* - this dependence is shown in Fig. HI For small values 
of a*, W re i takes large values and the extremization constraint is dominated by the (|, -)- 
mode causing the emitted fermions to be aligned with the rotation axis. As a* increases, W re \ 
decreases reaching the value one for approximately a* = 1 - now, both modes contribute equally 
and 9 max is pushed away from the 6 = tt value. For even larger values of a*, the (|, |)-mode 
starts dominating with the angle of maximum emission moving further away. 

In support of our argument, that the (|, |) and (|, §) modes predominantly determine 
the angular variation of the fermionic spectrum, in Tabled] we display the values of the energy 
emission rate at the angle of maximum emission as well as the value of the corresponding angle 
#max, for various values of the energy parameter Co>* and angular-momentum parameter a*. In 
each case, we display two values: the first one follows by taking into account the contribution 
of the two aforementioned modes and keeping terms only up to p = 3 in the sum of the angular 
eigenfunction (or, up to p = 4 for a* > 1); the second follows by keeping all terms up to p = 10 
and all partial modes up to I = 7/2. The values of the energy parameter have been chosen 
to lie in the low-energy regime and, at the same time, to display a non-trivial angular variation 
of the spectrum - it is worth mentioning that for all values smaller than = 0.5, the angle of 
maximum emission is constantly located at 9 = ir. On the other hand, the angular- momentum 
parameter a* scans a fairly broad range from a* = 0.5 to a* = 1.5. 

For the energy channel = 0.5, the agreement between the two sets of results is extremely 
good: the error in the value of the energy emission rate at its maximum reaches the magnitude 
of 3.5% at most, while the agreement in the value of 6 max is perfect. In agreement with the 
exact numerical results j30] where this energy channel was studied, for small values of a*, the 
emitted radiation remains very close to the rotation axis and only for values close to a* = 1.0 
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= 0.5 


= 0.6 


cu* = 0.7 


= 0.8 




approx. full 


approx. full 


approx. full 


approx. full 


a* = 0.50 


2.078 2.150 
-0.99 -0.99 


1.956 2.053 
-0.99 -0.99 


1.545 1.657 
-0.99 -0.89 


1.028 1.132 
-0.87 -0.74 


a* = 0.75 


2.316 2.406 
-0.99 -0.99 


2.054 2.168 
-0.82 -0.81 


1.693 1.827 
-0.66 -0.64 


1.232 1.402 
-0.57 -0.52 


a* = 1.00 


2.444 2.535 
-0.91 -0.91 


2.170 2.282 
-0.68 -0.67 


1.809 1.970 
-0.56 -0.54 


1.303 1.552 
-0.51 -0.46 


a* = 1.25 


2.488 2.567 
-0.81 -0.81 


2.090 2.205 
-0.61 -0.60 


1.590 1.784 
-0.52 -0.50 


1.019 1.338 
-0.46 -0.43 


a* = 1.50 


2.347 2.415 
-0.76 -0.76 


1.715 1.831 
-0.57 -0.56 


1.038 1.248 
-0.41 -0.42 


0.098 
-0.40 



Table 1: The approximated and full values of the energy emission rate (J7J) at the angle of 
maximum emission, in units of 10 _3 /r/ l , and the corresponding values of cos(# max ) for fermions. 



the emission starts showing a maximum at a gradually smaller angle. For = 0.6, the errors 
in the value of the emission rate and 6* max are at the level of 5% and 3% respectively, with the 
emission being peaked at an angle away from the horizon axis for a* > 0.75. For = 0.7, 
the error in the value of 9 maiX is still quite small ranging between 3% and 4%, whereas the 
error in the value of the emission rate at its maximum is now taking large values (7%-17%). 
Finally, for completeness, we show the energy channel of uo* = 0.8: although we have probably 
exceeded the range of validity of our approximation, the error in the value of 6 max remains less 
than 10%. 

The above comparison demonstrates that, for low values of the parameters and a* where 
our semi- analytic approximation is valid, the use of the two modes, the (|, |) and (§, §) ones, 
and the constraint ( 1351) can provide realistic results for the angular variation of the fermionic 
spectrum. This can consequently help to determine the value of the angular momentum of the 
black hole according to the proposal of [29| 130] . The results displayed in Table [1] confirm the 
behaviour found numerically for the energy channel = 0.5 [30], extend the set of values that 
could be used for comparison with experiment to additional low-energy values of Co>* and, finally, 
provide a very satisfactory semi-analytic approximation in terms of only two partial modes. 

4.2 Emission of Gauge Bosons 

Let us now address the emission of gauge bosons on the brane by the simply-rotating black 
hole. We will again focus on the low-energy regime as this is the energy channel at which 
the emission of gauge bosons is polarised along the rotation axis of the black hole. We will 
attempt to determine the main factors that contribute to this behaviour and, if possible, provide 
analytical arguments that justify it. 

4 The error in the value of 6 max is indeed quite small for all values of a* > 0.75. For a* = 0.5, we observe a 
significant deviation of (9 m ax from its actual value for the energy channels oj* — 0.7 and w* = 0.8. This is due 
to the fact that, for these specific values of the energy parameter and angular momentum , the mode (5, — 5) 
that we have ignored in our approximation is of the same order of magnitude as the (§, §) that we have taken 
into account. 
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COS^t7 max J 


p = 


0.017948 


-0.99 


0.029051 


-0.99 


0.067319 


-0.99 


p = 1 


0.017866 


-0.99 


0.028778 


-0.99 


0.066353 


-0.99 


p = 2 


0.017866 


-0.99 


0.028779 


-0.99 


0.066357 


-0.99 


p = 3 


0.017866 


-0.99 


0.028779 


-0.99 


0.066357 


-0.99 



Table 2: The differential energy emission rates at the maximum angle of emission and the 
corresponding angle for gauge bosons, for = 0.3 and n = 2, in terms of the angular- 
momentum number I and sum index p, and for three indicative values of a*. 

Following a similar strategy as in the case of fermions, we first investigate whether the 
infinite sum over the partial modes, characterised by (l,m), in Eq. (130]) can be truncated. 
By gradually increasing the value of / (and summing over all corresponding values of m), we 
looked for that value beyond which any increase in / makes no difference to the value of the 
energy emission rate at its maximum and of the corresponding angle. It turns out that, at the 
low-energy regime, this value is reached very quickly - this behaviour is clearly displayed by 
the entries of Table [2j In the upper part of the Table, we present the energy emission rate ([7]) 
at its maximum and the corresponding angle as we increase / from 1 to 3 and vary a* from 0.5 
to 1.5 in a random low-energy channel (a;* = 0.3). We observe that the value of the angle of 
maximum emission for positive-helicity (h = 1) gauge bosons is indeed 8 = ir, i.e. anti-parallel 
to the angular- momentum vector of the black hole, and that this value is not affected at all by 
adding any partial modes beyond the ones with 1 = 1. The energy emission rate also varies 
very little: its value at the angle of maximum emission is already reached for / = 2 and the 
difference from its value when only the / = 1 modes are taken into account is of the order of 
0.1% independently of the value of the angular- momentum of the black hole. We may thus 
conclude that the angular profile of the emission of gauge bosons at the low-energy regime is 
determined almost exclusively by the lower I = 1 modes: the sum over Z, therefore, in Eq. fT30|) 
can be replaced by the contribution of only its first term. 

We performed a similar analysis regarding the value of p in the sum in the expression 
for the angular eigenfunction, and we have found similar results displayed in the lower part of 
Table [2J The value of the angle of maximum emission is again not affected as terms beyond 
the first one (p = 0) are added. The actual value of the energy emission rate at the maximum 
angle is also very loosely dependent on p: as p goes from 1 to 2, the difference is of the order of 
10 _3 %, while the difference between the cases with p = and p = 1 is again very small, of the 
order of 0.5%. While, according to the above, the sum over p can be clearly truncated even at 
p = 0, to increase the validity of the subsequent analysis, we will also keep terms with p = 1, 
and thus write the analytic expression (|2"5"|) of the angular eigenfunction as 

iS^ar) = e°" x (1 + x) k ~ (1 - x) k+ [a + a x (1 + x)\ . (36) 

A final point that needs to be addressed is the contribution of the different m-modes. For 
/ = 1, we have three modes with m = +1, 0, —1 that have, nevertheless, a different weight factor 
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and thus a different contribution to the constraint (!30|) . A numerical evaluation of the weight 
factor (1311) . with Ta given in Eq. f ill I) , for these three modes, in conjunction with the value 
of the angular eigenfunction in each case, reveals that the contribution of the m = 1 mode to 
the constraint (1301) is almost two orders of magnitude larger than the one of the m = mode, 
and that in turn is larger by two orders of magnitude than the contribution of the m — — 1 
mode. Therefore, it is the I = m = 1 mode that effectively determines the angular profile of 
the emitted radiation. 

Then, the constraint f )30|) can take a particularly simple form. For I — m — 1 and h — 1, 
we obtain = and k + = 1, which then leads to the condition 

a (au - - — — ] + a 1 (1 + x) (au + — — — ] = . (37) 

\ 1 — X J \ 1 + X 1 — X J 

The above can be written as a quadratic polynomial in x, with solutions 

+ 4aj 




-ex = -(^ + ^)±Wl + 77^ + ? + 5- (38) 



If the above values correspond to extremal points in the regime x G (—1, 1), then they should 
satisfy the inequality |x ex | < 1. This in turn imposes constraints on the coefficients a and a\. 
As in the case of fermions, these coefficients, for a given set of numbers (h, I, m), are given solely 
in terms of the parameter aw. In Appendix \A.2\ we present the main steps for the derivation 
of the relations between the sum coefficients a p in the case of gauge bosons. There, it is found 



that, for the mode h = I = m = 1, 



(!) n 

a) au .„ 9au 



^- = -_ (3 + —). (39) 
W 2 v 20 



«o 

We substitute the above ratio into Eq. fl38l) . and demand that — 1 < x cx < 1. While the left- 
hand-side inequality is automatically satisfied for all values of au, the right-hand-side translates 
to 9(chj) 2 + 60au — 20 > that leads to the constraint au > 0.32. 

Therefore, for au < 0.32 no extremal points for the differential energy emission rate exist 
in the range x G (—1, 1). This quantity is thus monotonic and has global extremal points at 
the end points x = — 1 and x = +1. Substituting = and k + = 1 in Eq. ( 1361) . it is 
easy to see that, for x = +1, the angular eigenfunction vanishes, while, for x = —1, it takes 
its maximum value 2aoe~ au ' . As a result, the positive-helicity component of the gauge field 
is perfectly aligned in an anti-parallel direction to the angular-momentum vector of the black 
hole {6 = 7r), in agreement with the exact numerical results [30J. If au exceeds the value 0.32, 
a local maximum develops at an internal point of the range (—1,1), however, this remains 
subdominant to the global maximum at x — — 1 up to the value au ~ 0.85. Therefore, if we 
fix the energy channel to u* = 0.5, the maximum of the emitted radiation in the form of gauge 
fields remains aligned in an antiparallel direction to the angular-momentum of the black hole 
for all values of a* up to 1.7, in agreement again with the exact numerical results [30] . 



4.3 Emission of Scalars 

We finally address the case of the emission of scalar fields on the brane by a simply-rotating 
black hole. Although no useful information regarding the angular momentum of the black hole 
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Table 3: The differential energy emission rates at the maximum angle of emission and the 
corresponding angle, for = 0.3, n = 2 and a* = 1.5, in terms of the angular-momentum 
number I and sum index p for scalar fields. 

can be derived in this case, for completeness, we briefly discuss the main characteristics of the 
angular pattern of the scalar emission and the main contributing factors. 

In order to investigate whether it is possible again to truncate the sums over I and p, that 
appear in the constraint ( 1301) . we construct Table "21 The left-hand-side of the table displays 
the energy emission rate at the angle of maximum emission and the corresponding angle in 
terms of the angular-momentum number I. The energy channel = 0.3 has been chosen as an 
indicative case, the number of extra dimensions has been again fixed to n — 2, and the angular- 
momentum parameter is taken to be a* = 1.5 - this is the highest value of a* considered in this 
analysis, and the one for which the convergence of the sums over I and p is the most difficult 
to achieve. We observe that all modes beyond I = 2 add a contribution of order 0.01%, and 
thus can be safely ignored. But the difference between the values of the emission rate when 
all modes up to I = 1 and / = 2 have been, respectively, summed up is also very small, of the 
order of 1%. The value of 6* max has also been stabilised to n/2 when l max = 1. Therefore, in the 
context of our semi-analytic approach, the sum over I can be indeed truncated at / = 1 with 
no significant error. 

On the right-hand-side of Table EJ we keep all partial modes up to I = 2 for extra accuracy, 
and examine the convergence of the sum over p. The change in the value of the energy emission 
rate at the angle of maximum emission between the cases with p = 4 and p = 5 is of the order 
of 0.3%, while all higher contributions are an order of magnitude smaller. The value of (9 max 
has also taken the exact value of tt/2, therefore this sum can be safely truncated at p — 4. 

An exhaustive analysis of the values of the weight factors of the contributing partial 
modes (I, m) = {(0,0), (1, —1), (1,0), (1, 1)} for a variety of energy channels, cu* € (0.2 — 0.8) 
and angular-momentum of the black hole, a* G (0.5 — 1.5), reveals that the two most dominant 
modes are the (0, 0) and (1, 1) with the contributions of the other two being always two orders 
of magnitude smaller. Therefore, combining all the above results, the extremization constraint 
(13*01) for the case of scalar fields, takes the simplified form: 



Wrel E 4° 0) (1 + Xf 4° 0) (1 + ^ + Y+ 

4 4 / 

n—n r,—n V 



+ X 2(1 -x) 



. (40) 



p=0 g=0 

In the above, we have used that for the (0, 0)-mode, k- = k + = 0, while for the (1, l)-mode, 
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k- = k + = 1/2. Also, in this case, the relative weight factor is defined as W re i = Wqq/W\\. 
The expressions of the sum coefficients a p 00 ^ and a p 11 ^ for the two modes can be found at the 
Appendix IA.3I 

Let us consider individually the two dominant modes. Starting from the I — m — mode, 
we write its extremization constraint as 

3 

^2 \ auj a P + (P + 1) Vh ] (! + x f + auJ a 4(! + x) 4 = . (41) 

p=0 

This is a polynomial of fourth degree that in principle has four roots and, therefore, four 
potential extremal points. However, if we use the expressions of the a p coefficients for the 
/ = m = mode listed in Appendix IA.3[ we find that two of these roots are complex conjugates 
and one lies outside the range [—1,1]. Thus, the angular wavefunction of the I = m = mode 
has only one extremal point with respect to x — cos 9. This extremum is a minimum located 
at x = (9 = 7r/2) for small values of au that moves to positive values of x as au increases. 
However, the latter effect is actually an artifact of the truncation of the sum in the expression 
of the angular eigenfunction at a finite value of p. Even in our approximation where terms up 
to p = 4 are kept, we may see that the constant term of the polynomial (j41|) is given by a 
particular combination of the a p coefficients that due to multiple cancellations quickly tends to 
zero, namely 

[au a p + (p + 1) a p+1 } + au a 4 ~ + 0(auf . (42) 

p=0 

Had we kept all terms in the series expansions of the angular eigenvalue and eigenfunction, 
every subsequent term in the sum of fj42|) would cancel part of the remain of all previous ones 
all the way to infinity, thus ensuring that the x = is always an extremum of the I = m = 
mode. A simple numerical analysis then shows that this local extremum is the only one in the 
range (—1,1) and corresponds to a minimum. Due to the fact that k + = k- = 0, the I = m = 
mode reaches the same maximum value at the boundary points x — ±1. 

Moving to the next dominant mode I — m — 1, its extremization constraint reads 

4 

Y a P (! + X Y [ auj (1 - x 2 ) -x{p+l)+p] = . (43) 

p=0 

This is a polynomial of sixth degree whose six roots are potential extremal points. Substituting 
the a p coefficients for this mode from Appendix I A. 31 and performing a simple numerical analysis, 
one may see that the four roots are two pairs of complex conjugate numbers and one lies outside 
the range [—1, 1] leaving again only one root that may indeed correspond to a local extremal 
point of the angular eigenfunction of the I — m = 1 scalar mode with respect to x — cos 9. 
As in the case of the I = m = mode, the extremum is located at x = and moves towards 
positive values of x as the parameter au increases. We have again confirmed that the constant 
term of the above polynomial tends again to zero very quickly, i.e. 

+ p)a p ^ + 0{auf , (44) 

p=0 

signalling the fact that the x = is always an extremum of the angular eigenfunction of the 
I = m = 1 mode. The difference from the case of the / = m = mode lies in the fact that now 
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this extremum is a global maximum instead of a minimum with the angular eigenfunction of 
the I — m = 1 mode vanishing at the boundary points x = ±1 since k + = fc_ = 1/2. Let us 
briefly add here that a similar analysis of the remaining two scalar modes, I = l,m = and 
/ = — m = 1, shows that these follow the behaviour of the I = m = and I — m — 1 modes, 
respectively 

The exact numerical analysis of the emission of scalar fields on the brane by a simply- 
rotating higher- dimensional black hole [H] has revealed that the corresponding spectrum shows 
no angular variation for low values of the energy parameter cu* and of the angular-momentum 
number a*. Clearly, for a* = 0, the constraint ( 14"TI) is trivially satisfied and the I = m = 
mode shows no extremal points - note that, for the mode I = m = 1, the constraint ( |43l) 
still leads to a maximum at x = even at a* = 0. For low values of a;*, a careful analysis 
reveals that it is the I = m = mode that dominates over the others, therefore, for low a*, the 
spectrum remains spherically-symmetric. As a* starts increasing, the / = m = also develops 
an extremum at x — - it turns out that there is always a low-energy regime where the minima 
of the I = m = and / = l,m = modes exactly cancel the maxima of the I = m = 1 and 
I = — m = 1 modes leading again to a spherically symmetric spectrum, however, this energy 
regime becomes gradually more narrow. If we allow the energy parameter cu* to increase, too, 
then fairly quickly the I — m — 1 mode starts dominating causing the spectrum to exhibit 
maximum emission at x = 0, i.e. on the equatorial plane (9 = tt/2), in agreement with the 
exact numerical results [.14]. 

5 Discussion and conclusions 

One of the most exciting prospects of the theories predicting the existence of additional spacelike 
dimensions in nature and a low fundamental scale for gravity is the potential creation of higher- 
dimensional black holes from the collision of ordinary brane-localised particles. If the scale is 
low enough, the creation could in principle take place at ground-based particle accelerators and 
possibly be observed in the near future. The main observable signal is considered to be the 
emission of Hawking radiation on the brane in the form of ordinary Standard Model particles. 

During the study of the spherically-symmetric Schwarzschild phase, it was found that 
the radiation spectra of higher- dimensional black holes - even if we focus on the part of the 
emission that takes place on the brane where ourselves, the observers, are located - show a 
strong dependence on the number of additional spacelike dimensions that exist transversely to 
our brane. Therefore, the expectation was formed that the detection of the Hawking radiation 
spectra could lead to the determination of the number of extra spacelike dimensions in nature. 
However, if the angular-momentum of the black hole is taken into account - which generically 
is non-zero and seems to dominate almost all of the life of the black hole - this dependence on 
n is entangled with the dependence on the angular-momentum parameter a. 

During the spin-down phase of the life of the produced black hole, the emission exhibits, 
among other features, a strong angular variation in the radiation spectra with respect to the 
rotation axis of the black hole. It has been suggested (291 US] that this angular variation is 
the observable that could disentangle the dependence of the radiation spectra on n and a as 
it depends strongly on the latter while being (almost) insensitive to the former. It was found 
that, in the low-energy channel, the emitted gauge bosons become aligned to the rotation axis 
of the produced black hole while fermions form an angle with the rotation axis whose exact 
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value depends on the angular-momentum of the black hole. 

Attacking the problem of the angular variation of the Hawking radiation spectra in an 
exact way, and for all values of the parameters of the theory, is extremely challenging. It 
demands the numerical determination of both the radial and angular eigenfunction of the 
emitted fields as well as the numerical calculation of the angular eigenvalue that connects the 
corresponding equations. In addition, the angular pattern of the emitted spectra is formed from 
the contribution of an, in principle, infinite number of partial modes, numbered by the pair of 
angular-momentum numbers (l,m), each entering in the expression of the emission rate with 
its own weight (thermal and greybody) factor. Therefore, the use of the formal extremization 
constraint, that we have derived and which should determine the angles of maximum emission 
of all species of particles, seems rather unrealistic. 

Nevertheless, as the exact numerical analyses [29, 30J have shown, all the valuable infor- 
mation that we should deduce from the angular spectra are restricted in the low-energy regime. 
In this regime, one may use approximate techniques to solve the radial equation and thus 
determine the weight-factor of each contributing partial mode [TBI [19] . In addition, analytic 
formulae for the angular eigenvalue and eigenfunction exist [2SJ ESI EH ESI ESI BE] that allow 
us to study the problem of the angular variation of the spectra without resorting to complex 
numerical techniques. Combining the above tools in a constructive but critical manner, we 
were able to study the angular variation of the Hawking radiation spectra of fermions, gauge 
bosons and scalar fields in a semi-analytic way. 

Starting from the case of fermionic fields, the use of the analytic form of the greybody 
factor allowed us to compute the weight factor of each contributing partial mode. This, com- 
bined with a power series form for both the angular eigenvalue and eigenfunction, led to the 
isolation of the partial modes that predominantly determine the angular pattern of the cor- 
responding radiation spectra in the low-energy regime. Also, by demanding that the errors 
associated to the elimination of all higher-order terms were small, we were able to truncate 
the infinite sums in both the expressions of the angular eigenvalue and eigenfunction. At the 
end, we demonstrated that the contribution of only two partial modes, the (|, |) and (§,§), 
was more than adequate to provide approximate results for the value of the angle of maximum 
emission and of the corresponding emission rate that were within a range of 5% accuracy of the 
full results. Our study was completed by the derivation of the values of the above quantities, 
both in an exact and approximate way, for a variety of values of the energy parameter of 
the emitted fermionic field and angular-momentum parameter a* of the black hole, that could 
in principle be used for the determination of the angular momentum upon the observation of 
such a radiation spectrum. 

Whereas the angular variation of the radiation spectra of the emitted fermions is very 
sensitive to the value of the angular-momentum of the black hole - the larger the a parameter 
is, the larger the value of 6* max - the orientation of the gauge bosons emitted in the low-energy 
regime was found, by the exact numerical analyses, to be constantly aligned to the rotation 
axis of the black hole. Thus although it seems that no further information can be deduced from 
the study of the gauge bosons, we nevertheless performed the same analysis in an attempt to 
justify analytically the predicted behaviour. We demonstrated, by using a similar strategy as 
in the case of fermions, that a single mode, the I — m — 1, mainly determines the angular 
profile of the emitted gauge bosons. As its angular eigenfunction exhibits no extremal points 
up to au = 0.32, and even then these extrema remain subdominant to the global maximum 
at 9 = 0, 7r (for helicities h = =pl, respectively) up to au = 0.85, it is thus confirmed that 



the emission of gauge bosons in the low-energy regime will remain aligned to the rotation axis 
of the black hole for a wide range of the angular-momentum parameter. For example, gauge 
bosons emitted in the energy channel cu* = 0.5 will remain mostly parallel or antiparallel to the 
rotation axis up to the value of a* = 1.7, however, they will start deviating significantly from 
this behaviour for values of the angular momentum parameter of the black hole larger than 
this. 

For completeness, we have finally studied the case of the emission of scalar fields on the 
brane by a simply-rotating black hole. In this case, the exact numerical analyses have shown 
that the emission remains spherically-symmetric for low values of u and a and then, as either 
of the two parameters increases, the emission starts concentrating on the equatorial plane. One 
could thus assume that by looking at the emission of scalar fields, the equatorial plane, and thus 
the rotation axis of the black hole, could again be determined. Our analysis has shown that the 
angular profile of the radiation spectra of scalars in the low-energy regime is mainly determined 
by two modes, the I — m — and I — m — 1, the first having a minimum at 9 = n/2 and the 
second a maximum at the same point. For small values of a;* and a*, we have confirmed that the 
combination of these two modes creates indeed a "spherically-symmetric zone" in the emission 
where the two extrema exactly cancel each other. As either cu* and a* increases further, it is 
the I — m — 1 mode that starts dominating pushing the bulk of the emitted scalars towards 
the equatorial plane. Nevertheless, this transition becomes gradually and is finally realised for 
values of the parameters of the theory where our approximate techniques are not valid any 
more. 

Closing, let us note that our analysis in this work was based on the study of aspects, such 
as the existence of extremal points and relative magnitudes, of the spin-weighted spheroidal 
harmonics. These functions arise in a variety of problems, both in four dimensions as well as in 
the context of brane models, whenever the study of spin-s fields in a 4-dimensional spacetime 
with one angular-momentum component is performed. We thus envisage that the properties of 
the angular eigenf unctions revealed in this analysis as well as the analytic expressions of the 
angular eigenvalues and a p coefficients for fermions, gauge bosons and scalar fields presented in 
the Appendix will be of use in a variety of problems. Nevertheless, let us stress that the aspects 
of the particular problem studied here, i.e. the angular profile of the emission of Standard 
Model particles on the brane by a higher-dimensional black hole, could not be performed only 
by means of 4-dimensional tools: the particles emitted propagate on a brane embedded in a 
higher- dimensional spacetime, and this is reflected in the expressions of the greybody factors 
that determine to a great extent the weight factors of the individual partial modes. It is 
therefore the combination of both traditional 4-dimensional and brane techniques that has 
allowed us to analytically reproduce the angular distribution of energy emission, and hopefully 
provide the means for the determination of the angular momentum and axis of rotation of the 
produced black hole. 
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A Expansion Coefficients for the Spin- Weighted Spheroidal 
Harmonics 

The three-term recursion relations (|26|) - (|27|) can be used to determine the coefficients a p that 
appear in the expansion of the spin- weighted spheroidal harmonics hS\. For instance, for the 
first five expansion coefficients, we obtain 

A) Pi 7l f PIPO / yj f\ 

di = a , a 2 = ai a = a , (45) 

cto a.\ a.\ V^i^o Oil/ 

02 72 ( fcPiPo . P2J1 . 72/3o\ (AR s 

a 3 = a 2 ai = 1 1 a , (46) 

«2 Oi2 V a^ai^o a^cti ol^ol^ / 

/?3 73 ( fcfafkPo M2P0 IzfiiPo , 737i \ , A ^ 

a 4 = a 3 a 2 = 1 a . (47) 

«3 0J3 ya^o^i^o ol^ol-iOli c^a^ao a^ai^o a 3 «i / 

The (a p , (3 P , j p ) coefficients are given by Eqs. (|28|) and must be evaluated for each specific 
partial mode. 



A.l Fermions 

For the needs of our analysis, we determine the above coefficients for the fermionic modes 
(~, |) and (§,§)■ First, for the (~, ~)-mode, Eq. ( )2T|) leads to the following result for the 
corresponding eigenvalue (h = 1/2) 

CLOJ 1 1 . . o , 

Jn = H , 48 

2 22 3 27 

where we have used that, for this mode, k_ = and fc + = 1/2. Then, the (ct p , j p ) coefficients 
take the form 



2 1 fi 

< 2 ) = -2(p + l) 2 , / 3(V2) =p(p + 2) _4aa;( J 9+-)--(aa;) 2 , 7 < l /» = 2 au(p + 1) . (49) 



Then, the recursion relations (|26l)-(l27l) lead to the following relations between the first four sum 
coefficients 

-aw 1 + — , -^ = (aa;) 2 l + —— +..., (50) 



(1/2) 3 V 9 ' a (1/2) 



4 V2) 392 3 /. 1ST 



aa; ) I 1 + 777 aa; )+•••■ (51) 



(1/2) 729 ^ 'I 441 

The superscript (1/2) denotes that the above expressions hold for the case of the (|, |) partial 
mode. We have also given only the relations between the first four expansion coefficients since, 
as shown in section I4.1[ in the case of fermions, it suffices to consider only terms up to p = 3 
in the sum of Eq. ff25l . 
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For the mode (|, |), using the fact that now = 1/2 and fc+ = 1, we arrive at the 
following result for the angular eigenvalue 

au 29 . n9 . 

1^33=3 (au) 2 . 52 

222 5 125 v ; ( J 

In turn, the (a p , (3 p , j p ) coefficients take the form 

6 Q6 

a&ft) = - 2 (p + l)(p + 2) , /3f 2 ) = p(p + 4) - 4aw(p + | ) - ^ (aw) 2 , 7 f 2 ) = 2aw(p + 2) . 

(53) 

Finally, the above result into the following relations between the first four sum coefficients 

af /2) 6 / 4aw\ af 2) 4. , 2 / 34aw\ 
■j— = — au 1 , -4?i= r = -(aw) z H +•••, (54) 

. , 3 / 1588 \ 

M ( 4475 GW j + "" ( } 



A. 2 Gauge Bosons 

We now turn to the case of gauge bosons, and more particularly to the dominant mode with 
h = I = m = 1. Employing Eq. (122)) and (1281 . we find the following results for the angular 
eigenvalue 

11 
20 

and the (a p , f3 p , j p ) coefficients 

9(auj) 2 

a p = -2(p+l) 2 , (3 p = p(p + 3)-au(Ap + 3)-^^, lv = 2au{p + 2) . (57) 



iAn = —au — — (au) 2 , (56) 



We may then compute the relations between the different sum coefficients a p - although for our 
analysis we need only the relation between a\ and a Q , for completeness, we display again the 
relations between the first four coefficients 

a^ au 9au s a? , ,^ / 6 9au 



7i)=-T< 3 + ^- U + ^) + "" (58) 



(aw) 3 /65 1247 



a Q 



W 32 V 3 200 



1 aw + ... . (59) 



A. 3 Scalar Fields 

We finally study the case of scalar fields (h = 0). For our analysis, we will need the sum 
coefficients a p for the modes I = m = and I — m — 1. We start with the case with I = m = 0: 
employing again Eq. (124")) and (1281 . we find the following results for the angular eigenvalue 



(au) 2 , , 

oA 00 = , (60) 
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and the (a p , (3 p , j p ) coefficients 

a p = -2(p+l) 2 , f3 p = p{p + 1) - aw(4p + 2) - , 7p = 2awp, (61) 



that, in turn, lead to the following relations between the first five sum coefficients a 

(00) (00) / \2 / I \2 

a\ aw, a), (aw) /„ (aw) 



aw(l + — ), *=±-J- 2 + aw + ^-f , (62) 



a (oo) v 3 - fl (oo) 3 ^ 12 

4° 0) (aw) 3 / 65aaA ai° 0) _ (aw) 4 /115 29aw\ 

a (oo) 3 108 a (oo) 108 \8 + 3 J + ' [M) 

For the mode / = m = 1, a similar analysis leads to the following results for the eigenvalue 

(aw) 2 , , 

A n = 2 - ^-L , 64 
5 

the (a p , j p ) coefficients 

a p = -2(p + l)(p + 2) , /3 p = p(p + 3) - 4aw(p + 1) - , 7p = 2aw(p + 1) , (65) 

and the first five sum coefficients a p 

(H) (11) / \2 / / \2\ 

' aw. a, ; aw r / awr\ . s 

aw (1 + — ) , -2-r = 3 + aw + K — '- , (66) 



aj n) V 5 ' ' a W 5 V 15 

ai n) (aw) 3 / 103aw\ aj n) (aw) 4 /571 ^_ 43aw\ 
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